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Abstract 

Oh ' We analyze the rate of convergence towards self-similarity for the subcritical Keller- 

Segel system in the radially symmetric two-dimensional case and in the corresponding 
one-dimensional case for logarithmic interaction. We measure convergence in Wasser- 
stein distance. The rate of convergence towards self-similarity does not degenerate as 
we approach the critical case. As a byproduct, we obtain a proof of the logarithmic 
Hardy-Littlewood-Sobolev inequality in the one dimensional and radially symmetric two 
dimensional case based on optimal transport arguments. In addition we prove that the 
one-dimensional equation is a contraction with respect to Fourier distance in the subcrit- 
£> ' ical case. 
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We will concentrate on seeking decay rates towards equilibria or self-similarity profiles for 



o 

aggregation equations with linear diffusion in the fair competition regime. These models 
describe the evolution of a population of individuals which are diffusing by standard Brownian 
motion and attracting each other by a pairwise symmetric potential W(x). We focus on a 
logarithmic interaction potential W(x) = 2% log \x\, with x > 0. The Fokker-Planck equation 
governing the evolution of the probability density function p(t, x) associated to this particle 



system reads as 



at 



—Vp + 2 X p(Vlog\x\*p) 



t>o, xew\ (l.i) 



Due to translational invariance and mass conservation, in the rest of this work we restrict to 
zero center of mass probability densities, 

p(t, x)>0, / p(t, x) dx = 1 , / xp(t,x) dx = , 
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By fair competition, we mean that the dynamics of (|1.1|) are driven by a simple dichotomy as 
in the classical Keller-Segel system in two dimensions |24[ [30], [20], fl~2] , the modified Keller-Segel 
system in one dimension |151 [7] or the Keller-Segel model with suitable nonlinear diffusion 
in larger dimensions [9]. In all these examples there is a critical parameter which makes the 
distinction between global existence of solutions and finite-time blow-up. More precisely, we 
will discuss the modified one-dimensional Keller-Segel equation pT5| FT] : 

dtp = d 2 xxP + 2 x d x (pd x (log \x\ * p)) , t>0, iet, (1.2) 

and the radially symmetric two-dimensional classical Keller-Segel equation: 

d t (rp(t,r)) = ^d r (rd r p(t,r))+2 X d r [p(t,r)M[p(t)}(r)}, t>0, r € M+ , (1.3) 

where M[p] denotes the cumulated mass of p inside balls, 

M[p](r) = 2vr / p(s)sds. 
Jo 

Both equations (|1.2|) and (|1.3p exhibit a transition depending on the sensitivity coefficient x : 



• Subcritical Case.- For any < x < 1 solutions exist globally-in-time and they ap- 
proach a unique self-similar solution as t — > oo, see [201 EZ1 G3 [15] . 

• Critical Case.- For x = 1 solutions exist globally-in-time. There are infinitely many 
stationary solutions with infinite second moment. Solutions having finite initial second 
moment concentrate in infinite time towards the Dirac mass 5q \1U\ EJ [2S] . Solutions 
of infinite initial second moment close enough to a stationary solution converge to it as 
t->oo (Hj- 

• Supercritical Case.- For any x > 1 smooth fast-decaying solutions do not exist glob- 
ally in time [mi ei na us]. 

The critical parameter x = 1 can be obtained from two formal computations at this stage. 
The evolution of the second moment satisfies in both cases the relation: 



1 d f o . 

/ \x\ p(t,x) ax = 1 

2dtJ R l ' n ' 



X' 



This implies that for x > 1 solutions will necessarily blow-up before the second moment 
touches zero. On the other hand, the Keller-Segel equation (jl.ip is equipped with a free 
energy (entropy minus potential energy), 

F ^ = 1TtI P( x ) l °gP( x )dx + x p(x)log(\x-y\)p(y)dxdy. (1.4) 

It is formally decreasing along the trajectories 

\ 2 

dx. (1.5) 



d 



di F[p(t)] = - J p(t,x) 



V ( — log p(t, x) + 2x log \x\ * pit, x) 



Moreover, it was shown in |2U[ [T2] that for % < 1 the free energy estimate from above implies 
an a priori bound in the entropy part of the functional which is at the basis of the construction 
of global-in-time solutions. This was achieved by using the Logarithmic-HLS inequality [Hll7| 
which relates the entropy and the interaction part of the functional. 

Nontrivial equilibrium profiles or critical profiles, only exist for the critical parameter 
X = 1. They are solutions to the following Euler-Lagrange equations: 

(i'(x) + 2n(x)8 x (log |x| * (i(x)) = , (1.6) 

^rp'(r) + 2p(r)M[p](r) = 0, (1.7) 

resp. in dimension N = 1 and in dimension N = 2 with radially symmetry. In fact, we have 
an explicit formulation of the stationary states, 

^ = ^nW' iV = 1 ' 2 - (L8) 

This coincides with the equality cases in the Logarithmic-HLS inequality. 

In the subcritical case x < lj solutions are known to converge to unique self-similar profiles 
|12| . For studying convergence towards self-similarity, it is generally useful to rescale the space 
and time variables in the subcritical regime x < 1- The Keller-Segel system rewrites as 

^ = ^Ap + 2 X V-[p(Vlog\x\*p)} + S/-[xp] , t>0, xeR N , (1.9) 

at JS 

and the free energy is complemented with a quadratic confinement potential: 

-^rescM = F[P] + I f \X\ 2 P&) <& • (1-10) 

2 Jr n 

Due to the change of variables, self-similar solutions correspond to equilibrium solutions of 
(jl.9p . The rate of convergence towards equilibrium for (II. 9p in the subcritical case was recently 
studied in |llj where the same rate as for the heat equation was obtained for small mass. 

Let us finally mention that both (11. ip and (jl.9p are gradient flows of the free energy 
functionals (|1.4p and (jl.lOp respectively, when the space of probability measures is endowed 
with the euclidean Wasserstein metric Wi- We refer to the seminal papers \12\ 129] and to [I] 
for a general theory. For instance, we can write (jl.ip in short as 

p(t) = -V W2 T[p{t)} . (1.11) 

This assertion was made rigorous in [7], where the variational minimizing movement scheme 
PQ was shown to converge for (jl.ip . This fact allows us to consider a way of measuring 
the distance towards equilibrium or self-similarity intimately related to the evolution due to 
(jl.lip . In fact, we will show that optimal transport tools are key techniques to describe this 
behavior at least in the one dimensional case fll -2f) and in the radial case in two dimensions 

do}. 

In order to investigate further the bounds of the free energy functional leading to the 
dichotomy discussed above and the characterization of the critical profiles, the Logarithmic- 
HLS inequality proved in [H |T7j is essential. In Section 2, we show an alternative proof 



based on optimal transport tools in the one dimensional case and in the radial case in two 
dimensions. There is another recent proof of this inequality with sharp constants in the two 
dimensional case by fast diffusion flows |16j . The Logarithmic-HLS inequality can be restated 
with our notations as: 

Theorem 1.1 (Logarithmic HLS inequality). Assume N = 1,2. The functional J 7 is bounded 
from below. The extremal functions are uniquely given by (jl.8p up to dilations in the set of 
probability densities with zero center of mass. 

In short, we demonstrate that any critical point of the free energy is in fact a global 
minimizer. This is a property which holds true for convex functionals, although the functional 
J- is not displacement convex in the sense of McCann |27j . 

The ideas behind the proof of the sharp Logarithmic-HLS inequality allow us to tackle the 
rate of convergence in W 2 by similar methods for the rescaled version (|1.9|) in one dimension 
and for radial densities in the two dimensional case provided x < 1. We prove in Section 3 
the following result. 

Theorem 1.2 (Long-time asymptotics). Assume that N = 1,2 being the initial data po 
radially symmetric if N = 2. In the subcritical case % < 1, solutions of (|1.9p in the rescaled 
variables converge exponentially fast towards the unique equilibrium configuration v. More 
precisely, the following estimate holds true 

jW 2 (p(t),uf <-2W 2 (p{t),uf . 

Surprisingly enough, the rate of convergence that we obtain does not depend on the 
parameter \- 0m - estimate is uniform as long as \ remains subcritical and is equal to the 
rate of convergence towards self-similarity for the heat equation. This is due to the fact that 
entropy and interaction contributions cancel each other, and only the confinement contribution 
remains yielding a uniform estimate. Although convergence is likely to be uniform, notice that 
the asymptotic profile becomes more and more singular as \ ~~* 1~ j as shown by the simple 
second moment identity 

\a\ 2 v(a)da = l-x- (1.12) 



Finally, we devote Section 4 to propose an alternative method of measuring the distance 
towards self-similarity in the one dimensional case. We make a connection between the one 
dimensional modified Keller-Segel model (jl.2p and certain Boltzmann-like equations used in 
granular gases and wealth-distribution models, see [181 [2T] and the references therein. This 
connection is due to the fact that (|1.2p can be written in Fourier variables like the referred 
Boltzmann equations. Following the ideas of [18] we prove that equation (]1.2p is indeed a 
contraction for the so-called Fourier distances defined in Section 4. 

Theorem 1.3. Assume x < 1 an d the initial data have finite second moments. The one- 
dimensional Keller-Segel system (II. 2p is a contraction for the distance d\. It is a uniformly 
strict contraction in the rescaled frame, with a contraction factor which does not depend on \- 



2 An alternative proof of the logarithmic HLS inequality 

2.1 Preliminaries on Optimal Transport Tools 



Let (j, and p be two density probabilities. According to [13, 26J there exists a convex function 
ip whose gradient pushes forward the measure p(a) da onto p(x) dx: Vip# (M°0 da) = p(x) dx. 
This convex function satisfies the Monge-Ampre equation in the weak sense, 

p(a) = p(Vip{a)) det D 2 ip(a) . 

Regularity of the transport map is a big issue in general. Here we will use the fact that the 
Hessian measure det HD 2 ip(a) can be decomposed in an absolute continuous part det aD 2 ijj{o) 
and a positive singular measure [31| Chapter 4]. In particular we have det HD 2 ip(a) > 
det AD 2 ip(a). The formula for the change of variables will be important when dealing with 
the entropy contribution. For any measurable function U, bounded below such that U(0) = 
we have [23 EI] 

17 W,))* - f RN U ( Hi j^ R ) de t , D V(«)<fa. (2.1) 

In fact this paper will only be concerned with the one-dimensional case, and the two- 
dimensional radial case. The complexity of Brenier's transport problem dramatically reduces 
in both cases. In dimension one, the transport map <p' is explicitely given by: ip'(a) = 
X o A" 1 (a) where X and A denote respectively the pseudo-inverse cumulative distribution 
function of the densities p and p. The singular part of the positive measure ip" corresponds 
to having holes in the support of the density p. 

In the two-dimensional radial case, the Brenier's map can be expressed as the one- 
dimensional transport between the densities 27rp,(a)a d ~ 1 da and 2irp(r)r d ~ l dr. The deter- 
minant of the Hessian is given by 

d«., B V(.) = ^(«f(»), 

where the derivative of (ip 1 ) 2 has to be understood in the distributional sense. 

The following Lemma will be used to estimate the interaction contribution in the free 
energy, and in the evolution of the Wasserstein distance. For notational convenience we 
denote the convex combination of a and b by [a, b]t = (1 — t)a + tb. 

Lemma 2.1. Let K : (0, oo) — > R be an increasing and concave function such that lim^^o F(z) - 
— oo. Then 

k(J ^"{[a,b] t )dt\ >^ K{^ c ([a,b] t )) dt. (2.2) 

Equality is achieved in ()2.2p if and only if the distributional derivative of the transport map 
ip" is a constant function. 

Analogously in the two-dimensional radially symmetric case we deduce 

I<( j det H D 2 i> ([a, b] t )dt) > I K (det A D 2 ip([a,b} t )) dt . (2.3) 

Equality is achieved in (|2.3[) if and only if ip' is a multiple of the identity. 



Proof. We have on the one hand ip" > tp" c . We next use the concavity of K to conclude. 
Equality occurs if t/j" is absolutely continuous and if ifj" c is constant. In the two-dimensional 
case we use det nD 2 ^(a) > det aD 2 iP (a). D 

Optimal transport is a powerful tool for reducing functional inequalities onto pointwise 
inequalities (e.g. matrix inequalities). We highlight for example the seminal paper by McCann 
|27| where the displacement convexity issue for some energy functional is reduced to the 
concavity of det ' . We also refer to the works of Barthe [21 [3] and Cordero-Erausquin et al. 
|19j . We require simple pointwise inequalities which are extensions of the classical Jensen's 
inequality. 

Lemma 2.2. We have the following convex-like inequality for some exponent 7 > and any 
positive u,v,a,j3, 

a (^iY.j^iY< {a+m (^±^i). 2 ,. ( , 4) 



Equality occurs if and only if u = v = 1. The continuous version reads as follows. For any 
measurable function u : (0, 1) — > (0, +00): 

a(f u(t)dt\ -P([ u(t)dt\ < (a + fi) f (u(i)) -7 dt - 2/3 , (2.5) 

Proof We only prove (12. 4j) . The continuous version (]2.5p is obtained by an approximation 
procedure. We introduce the auxiliary function J defined as follows. 

J(u,v) = (a + I3)( \-a\—-J + /3 

Clearly, J diverges towards +00 as u — > or v — > 0, and as u — > 00 or v — > 00, and J 
is bounded below. Then there exists at least one critical point. Any critical point (uq,vq) 
satisfies 

a + /3 - 7 -i , a ( uq + vq\~~ 1 ~ i _ /3 ( u Q + v x 



-^^^ +^2 ^ — 2 — ; + m-^-j ^ 

Q + /3 _ 7 _i a/n + w \ 7 1 /3 { u + v Y< l 



Hence «o = ^0 and 



^-- 1 +^- I +k I =»- 



2 u 2 u 2 

We conclude that u^ 7 = Uq. Therefore the unique critical point of J is (1, 1). □ 

2.2 The one-dimensional case 

The novelty here is contained in the proof of the logarithmic HLS inequality. This brings 
no information by itself since the uniqueness of the extremal functions is already known |17] . 
We show below that the logarithmic HLS inequality is a simple consequence of the Jensen's 

6 



inequality. However our proof relies on the existence of a critical point of the free energy T ' . In 
short, we demonstrate that any critical point of the free energy is in fact a global minimizer. 
This is a property which holds true for convex functionals. However the functional here is not 
convex. 

Our first Lemma is a reformulation of the Euler-Lagrange equation for the extremal func- 
tion (01) . 



Lemma 2.3 (Characterization of extremal functions). The critical profiles satisfy the follow- 
ing identity, 



Kp) 



/ / MP ~ ta )KP — tq + Q) dtdq . (2.6) 

Jr Jo 

In the subcritical regime x < 1, the equilibrium in the rescaled frame satisfies the following 
identity, 

"(P)= [ [ (x + 4-)v(p-tq)v(p-tq + q)dtdq. (2.7) 

JqeR JO \ z J 

Proof. The formulation (j2.6|) is equivalent to integrating once the equation for the critical 
profile. We integrate equation (11. 6ft against some test function tp. 



y'(p)n(p) dp = 2 I/ V>{x)n(y) dxdy 

xi x ~ y 

K x )Kv) dxdy 

K x-y 

l 

<f' ([x, y]t) (J,{x)n(y) dtdxdy 

xl JO 



II 

J Jr 



¥>'(p) 



/ / H(p - tq)^ip — tq + q) dtdq > dp , 
JrJo J 



where we have finally used the change of variables: (x,y) t— > (p = [x,y]t,q = y — x). This 
holds true for any derivative <p', so we obtain identity fj2.6[) up to a constant. Since both 
sides of (12.6|) have mass one, the constant is zero. The identity (12.7P is obtained in a similar 
way. D 

Proof of Theorem \l.l\ Applying the change of variables formula (12. ID for x = ip'(p), the 
functional J- rewrites as follows, 

F\p] - F\n\ = /log ( -^U n(a) da+ ff log \ip'{a) - ^(b)\fx(a)^b) dadb - F\jm] 

Jr \wLwJ JJrxr 

'log «») n{a) da + // log ( ^ (a) ~ f {b) ) m (o)m(6) dadb 

R J JRxR V a — J 

log (V4'c( a )) A* (a) da + // lo g ( / ^"([a> b ]t) dt J n{a)n{b) dadb 



xl VO 



Using Lemma 12.11 for K = \ogz which is increasing and concave, we deduce 

F\p] ~ Jfy] > - /log (C(p)) Kp) dp+ [[ [ log (C([«, b] t )) n(a)fi(b) dtdadb 



log(V&(p))/*(p)dp 

log «c(p)) { / M(P - *9)M(P - *9 + 9) ^idg Wp = . 

I LJqeEJO J 



+ 

Equality arises if and only if the transport map ip" is a constant function. Such a map 
corresponds exactly to the dilations of the critical profile \i. □ 

It is possible to extend Theorem 11.11 to the rescaled energy T vesc (|1.1U|) . 

Theorem 2.4 (Logarithmic HLS inequality with a quadratic confinement). Assume N = 1,2 
and x < 1. XTie functional «F res c * s bounded from below. The extremal functions are unique 
in the set of probability densities with zero center of mass. 

We give below the main lines of the proof following a direct argument analogous to the 
proof of Theorem 11.11 Note that the uniqueness of the extremal functions in dimension N = 1 
or in dimension N = 2 in the class of radially symmetric densities is a consequence of Theorem 
Ol 

Sketch of proof of Theorem \2.4\ The key point consists in replacing the Jensen's inequality 
with the following convex-like inequality. For any positive u, v, a, (3 the following inequality 
holds true. 

fu + v\ „fu-\-v\ , „. /log w + logtA 

a]OS \~2~) +P [~2~) ^( a + 2 ^[ 2 )+^ ( 2 - 8 ) 

Equality occurs if and only if u = v = 1. It reduces to the usual Jensen's inequality when 
/3 = 0. The proof of (|2.8p is analogous to Lemma [2.21 The proof of uniqueness for the extremal 
functions of .F resc is a mixture between the proofs of Theorem 11.11 and Theorem 11.21 □ 

2.3 The two-dimensional case 

We restrict to radially symmetric functions in the two-dimensional case due to decreasing 
rearrangement |25|, HI [T7]. We recall the Newton's theorem for Poisson potential: the field 
induced by a radially symmetric distribution of masses outside a given ball is equivalent to 
the field induced by a point at the center of the ball [25]. Equivalently it reads 

1 f 2jT 

- / log (r 2 + s 2 - 2rs cos(0)) dO = 2vr log max(r, s) . (2.9) 

2 Je=o 

As a consequence we can rewrite the functional T simpler under radial symmetry: 
TT^lfil = o / P( r ) l°g(/°( r )) rdr + x / p(r)M[p] (r) log(r) rdr . 
The following characterizations are direct consequences of (jl.7p and (jl.3p . 



Lemma 2.5 (Characterization of extremal functions under radial symmetry). The critical 
profiles satisfy the following identity 

1 f + °° 1 

-u(b) = 2 n(a)M[u](a)-da. (2.10) 

2 J b a 

In the subcritical regime % < 1, the radially-symmetric equilibrium satisfies the following 
identity 

-v(b)= u(a) (2xM[v](a)- + a) da. (2.11) 

We are now ready to examinate the logarithmic Hardy-Littlewood-Sobolev inequality in 
the two-dimensional radial setting. 

Proof of Theorem \l.l[ We apply the change of variables formula (12.ip for r = ip'(a) to get: 
^-F[P\ = \I M(o)kg ( , . M ffi,, A ada + 2 f /x(o)M[ M ] (a) log (V'(a)) ada, 

where we have used M[p](r) = M[/j,](a). We have consequently, 

1 1 1 /" 

^-^[p] - TrHlA = ~ o / ^( a ) lo S ( det a£> 2 ^(o)) a da 

2vr 27T 2 J R 



+ 1 fi(a)M[ f i}(a)\og( ^/^ Jo da. (2.12) 



The last contribution of (|2.12p can be evaluated using Lemma 12. II 

n(a)M[fj](a)log ( / (detH-D 2 ^(o)) —db) ada 

>/ / n(a) M[n](a) log (det A D 2 tfj(b)) — dbda 
Jr + Jo a 

= i log (det A D 2 i)(b)) 12 i n(a)M\p](a)- dalbdb. 

We obtain from the characterization (I2.10J) F[p\ > Tip]. Again equality occurs if and only if 
the transport map ip' is a multiple of the identity. □ 

2.4 Obstruction in dimension higher than three 

We explain in this Section why the above strategy fails to work in dimension higher than three, 
even in the radially-symmetric setting. A first remark is that Newton's Theorem is not valid, 
since the logarithm kernel is not the fundamental solution of the Poisson equation, although 
this is not essential as shown in dimension one. It turns out that our strategy works fine for 
any interaction kernel W(x) = \x\ k /k, for k € (—N,2 — N], The case k = corresponds to 
W = log \x\. The case k = —N is critical for integrability reasons. The case k = 2 — TV is 

9 



exactly the harmonic case for which the Newton's Theorem holds true. We refer to [13] for 
details in the case k € (-N, 2 — N]. Hence the case k = is out of range when N > 3. We 
sketch below where some obstruction appears when N = 3. 

The identity which generalizes (|2.9[) reads as follows. If r > s we have 



1 f n 1 /s\ 

- / log (r 2 + s 2 - 2rs cos(0)) sin(0) d9 = -H ( - ) + 2 log r , 

where -ff is defined as follows for t € (0, 1), 

aw= i(l+2! fcg(1+t) _>(lz5! l0g(1 _ _ 2 . 

To continue our strategy, it is required to decouple the variables r and s, and more precisely 
to make the quantity r — s appearing. As a matter of fact this is homogeneous to the 
determinant (under radial symmetry), which is the key quantity to look at in dimension 
higher than two. Therefore we seek a convex- like inequality 

^(t)>a + /31og(l-t iV ) , 

where a and /3 are suitable constants determined by zero and first order conditions. If we 
denote tp(t) = log(l — t ), this is equivalent to say that H o tp^ 1 is convex. However simple 
computations show that it is indeed a concave function. In the case of an interaction kernel 
having homogeneity k £ (-N, 2 — N] we show in [14] that the corresponding function H oip^ 1 
is convex. 

3 Exponential convergence towards the self-similar profile 

3.1 The one-dimensional case 

To illustrate the strategy of proof of Theorem 11.21 we show a formal computation in the 
critical case X = 1. Up to our knowledge, the regularity of solutions under very weak assump- 
tions is still an open problem In particular it is not known whether the solutions satisfy the 
identity (J1.5J1 or not. So the following computation is questionable because the velocity field 
d x (log p(t,x) + 21og|x[ *p(t,x)) is not clearly defined in L 2 (p(t,x)dx). 

We compute formally the evolution of the Wasserstein distance to one of the equilibria 
(jl.8p in the critical case x = 1- Notice that equilibria are infinitely far from each other with 
respect to the Wasserstein distance [8]. Using the gradient flow structure with respect to W2, 
one obtains the following formula for the derivative of F(t) = W2(p(t),p) 2 , see |3H Chapter 
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8] and [I]. 

1 d 
Ydt 



F(t) 



>'(t, x) — x) (d x (log p(t, x) + 2 log |»| * p(t, x))) p(t, x) dx 

<f/(t,x)-4/(t t v) 



< 



(j)"(t,x)p(t,x) dx + 
{%l)"{t,a))~ 1 p{a)da + 
(ip"(t 1 a))~ 1 p{a)da + 



x-y 



-p(t,x)p{t,y)dxdy 



^(t,a)-ip>(t,b) 
a — b 



p{a)p{b) dadb 



xl JO 



(ip"(t, [a, b] s )) p,{a)p,{b) dsdadb . 



We recognize the characterization f|2.6[) . Hence, we have at least formally F'(t) < 0. Observe 
that the Lemma 12.11 has been used with K(z) = — z _1 . 

The same strategy is valid in the subcritical case x < 1 f° r which we know that solutions 
are regular enough to ensure the validity of the computations. As a matter of fact, the density 
p(t, x) is everywhere positive and thus tp" is absolutely continuous. On the other hand the 
dissipation of energy is well-defined and the dissipation estimate (jl.5p holds true [12] . 



Proof of Theorem \1.1A We compute the evolution of F(t) = W2(p(t), v) 2 : 

2Jt F ^ = I (< ^ (f ' ^ ~ X) V x ( l0g P(t ' X) + 2X bg |X| * P(t ' X) + 2 ) ) / ° (t ' ^ dx 
= - f <j>"{t, x)p(t, x)dx + X If feL^M pfo x)p{ t, y) dxdy 

JR J JRxR X — y 

z iilxl 

+ 2 / <^'(t, x)xp(t, x) dx + 1 — x — / |x| 2 /)(£,x) dx , 
Jm. Jr 

where we have used the fact that the center of mass is zero to double the variables. We rewrite 
each contribution using the reverse transport map i\)': 



1 d „, s 

-—Fit) 
2dt K ' 



(i/)"(t, a)) v(a)da + x 



^{t,a)-^{t,b) 
a — b 



-i 



v{a)v{b) dadb 



,iP'(t,a)-iP'(t,b) 



\a — b\ — — — r v ' - u(a)u(b) dadb 

a — b 



< 



+ 1-X- / \ip'{t,a)\ 2 v(a)dx + 2 / aip f (t,a)u(t,a)da 
a-b\ 2 ^ rl 



X + 



{ip"(t,[a,b] a )) L ds-\a-b\< 



v{a)v(b) dadb 



(ip"(t,a)) u(a)da + 2 / \a\ 2 u(a) da — / \ijj'(t, a) — a\ 2 u(a) da, 
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where we have used that the second moment of the stationary state is explicitly given by 
(|1.12|) . Applying now (|2.5p for 7 = 1, a = x an d /3 = \a — b\ 2 /2, we deduce 

~-fF{t) < ~ / Wit, a) - a\ 2 u{a) da = -W 2 (p(t),v) 2 = -F(t) , 
z at j R 

giving the desired inequality. □ 

3.2 The two-dimensional radially-symmetric case 

Proving convergence towards a self-similar profile in the rescaled logarithmic case under radial 
symmetry goes as previously. 

Proof of Theorem M.SX The virial computation reads equivalently 

77-77 / p(t,r)r 3 dr = - / r I -d r log p(t,r) + 2xM[p](t,r)- +r ) p(t,r)rdr 
2dtJ R+ J R+ \2 - ■ r J 

p(t,r)r dr — 2\ I M[p](t,r)p(t,r)r dr — I p(t,r)r dr 



+ 

3 j„ — / ,,f„\„3 j„ I „/ + „\„3 



=1 — x — / p(t,r)r dr = v{a)a da — j p(t,r)r dr 

Jr + Jr + Jr + 

We compute again the evolution of the Wasserstein distance F(t) = W 2 {p{t), u) 2 . 
\j t F{t) = f (tf'(r)-r) ( ~d r log p(t,r) + 2 X M[p](t,r)^ + rj p(t,r)rdr 

=- r(f>'(r)d r p(t,r)dr + 2x (f>'(r)M[p}(t,r)p(t,r) dr - <f>' (r) p(t , r)r 2 dr 

2 Jr + Jr + Jr + 

+ 



+ / u(a)a 3 da- \ p(t,r)r 3 dr + 2 \ cf)'(r)p(t,r)r 2 dr 

-1/2 



<2 X f ( f det D 2 ip(b)^db\ M[v]{a)v(a)ada 

-/ L T^TTTv ) v{b)bdb- I ( I det D 2 4>{b)^db) a 2 u{a)ada 
Jr+ \det D 2 if)(b) ) Jr + \J a 2 J 

+ 2/ v{a)a 3 da-l \<f>'(r) - r\ 2 p(t,r)r dr . 

Jr+ Jr+ 

The last step in the inequality is a consequence of the arithmetic and geometric means in- 
equality: —d r (r<f/(r))/r = —cj)'(r)/r — 4>"(r) < — 2((j)" (r)cf)' (r) /r) 1 ' 2 . Next we use Lemma [2~! 
to handle the interaction contribution. More precisely, we choose 7 = 1/2, a = 2xM[v\{a 
and j3 = a 2 . One gets finally: 

~F(t)<-F(t)- f [det D 2 mY l,2 v(b)bdb 
* at Jr + 

{2 X M[u]{a) + a 2 ) (detD 2 ifj(b)y 1/2 ^v(a)adbda. 
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We conlude using characterization (|2.1ip . D 

4 Contraction in the one-dimensional case 

The aim of this Section is to point out the peculiar structure of the modified one-dimensional 
Keller-Segel system (|1.2p . 

Lemma 4.1. Equation (jl,2p rewrites in Fourier variables as: 

dtp(t, = |£| 2 (~p(t, O + xJ p(t, a£)p(t, (1 - a)0 da) . (4.1) 



Proof. We test equation (|1.2|) against exp(i£x): 

0_ 

Ot 



?(<! = / ( <W(*, a:) + 2xd x I p(t, x) I p.v. - j * p(t, x] 



e i£x _ p i£y 



= -\£\ 2 p(t,0-Xi£ // P(t,x)- — p(t,y)dxdy 

= ~ \Z\ 2 p(t, + Xl^l 2 If Pit, x) ( f 1 e«[**l- da) p(t, y) dxdy 

= -\i?P^i) + x\i? f] f p(t,x)e^-^ x dx) ( [ p(t,y)e^y d y) da, 

JO \JR J \JR J 

which gives the desired formulation. □ 

According to (|4,ip the information propagates from lower to higher frequencies. The 
evolution of p(t,£) requires the knowledge of lower frequencies |£'| < |^| due to the integral 
contribution. This is of particular importance for designing a numerical scheme. Indeed there 
is no loss of information after truncation of the frequency box. 

Remark 4.2 (Analogy with ID Boltzmann). It is worthy to metion that the integral operator 
in the right-hand- side of (14. ip is reminiscent of the homogeneous Boltzmann equations in ID 
used for granular gases fig]/ or wealth distribution models [21] in Fourier variables. 

Remark 4.3 (Evidence for blow-up in the supercritical case). We can directly notice the 
occurence of blow-up when x > 1 from (|4.ip . Observe that for |£| ^C 1, the right-hand-side is 
equivalent to: 

dtp(t, ~ I?! 2 {-P(t, 0) + X p(t, 0) 2 ) = |C| 2 (-1 + X) ■ (4-2) 

This argues in favor of blow-up at low modes although misleadingly. We have plotted in Figure 
U\ numerical simulation of (14. ip in the supercritical case. Observe that blow-up arises for 
|£| S> 1, on the contrary to the misleading heuristics (|4.2p . The integro- differential equation 
(14. 1|) makes perfect sense even in the supercritical regime X > 1 after the first blow-up event. 
However the outcoming fonction p(t, £) is no longer the Fourier transform of a probability 
measure. In fact the blow-up time coincides with the formation of the first dirac mass, namely 
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Figure 1: Simulation of (|4.ip in the supercritical case x > 1 f° r successive times. The blow-up 
time is plotted in bolded dark. 



when the frequency distribution p(t, £) is flat at infinity. This contradictory intuition is similar 
to the proof of blow-up based on the virial identity: the second momentum provides information 
at infinity but is used to prove blow-up which is a local behaviour. 

Recall the definition of Fourier distances [T5] as they have been introduced for the analysis 
of the Boltzmann equation. 

Definition 4.4 (Fourier distances). Let p\, p 2 being two probability measures having the same 
center of mass. The d\ — distance is defined as follows: 



d l ( Pl ,p 2 ) = snp{\C\- 1 \pi(0-p2(0\} 



(4.3) 



Proof of Theorem \1.3l First notice that supremum in (J4.3I) is attained in R\ {0}. Clearly we 
have |/3i(£) - /02(C) I < 2 an d 

Pi(0 - /82(e) ~ (J M 2 M*) " P2(x)\ dx\ \i\ 2 /2 as i -> 0. 

We denote F{t) = di(p-i{t) , p 2 {t)) and h(t,£) = |£| -1 (pi(*,£) - hi^O)- We multiply the 
difference between the two equations (|4.ip by sign(/i(£,£)), 



^IM<.OI = ie| 2 (-W*,OI+Xsign(pi(*,0-/32(t,e))^(*,0) 



where 



A(t,0 = \Z\- 1 I pi(t,a0pi(t,a-v)0do--\£\- 1 [ p 2 (t,a0f>2(t,(l-o-)0do-. 
Jo Jo 
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The self-attraction contributions are handled as follows \18\ Th. 6.3]: 

\A(t,0\ <|^ / \pi{t,crS)-P2(t,<rO\\pi(t,{l-<r)S)\d<r 

Jo 

+ ler 1 / \pi(t,(l-a)0-P2(t,(l-a)0\\p2(t,aO\da 
Jo 



< di ( Pl (t),p 2 (t)) I (a + (1 - a)) da = F(t) . 
Jo 



dt\h(t,0\<\Z\ 2 (-\h(t,d)\+xF(t)) 



\h(t + e,0\<e-W \h(t,0\ +X 1-e-^l sup F(t + s), 

V J s€(0,e) 



We obtain finally 
We deduce 



\h(t + 6,01 " F(t) < (1 - e - e| « |2 ) -F(i) + X sup F(t + s) , 
V y V se ^ / 

li msU p F(t + e) ~ F( ^ < (x " 1) ( liminf |C*(* + e)| 2 ^ ^(t) , 
e^o+ e V e->0+ / 

where |£*(i)| denotes the lowest frequency moduli for which the supremum is attained in 
F(t) = sup \h(t, £)|. We have used the continuity of F to pass to the limit. Therefore we get 
a contraction estimate as soon as x < 1- There is no explicit rate since we do not know how 
to control |£*(t)| from below. 

We also obtain a uniform strict contractivity in self-similar variables. The Keller-Segel 
equation (II. 9p writes as follows in Fourier variables: 



d tP (t, = |£| 2 (-p(t, Z) + xf P(t, aZW, (1 " <r)Z) da) - ^p(t, £) . 
We proceed as above to get: 

9t \h(t,0\ =ld 2 (- \h(t,0\ + xsign(pi(i,6 - fcftOMCU)) 
-^(\h(t,0\)-\Kt,0\- 
We integrate along characteristics and argue as previously, 

\h(t + e,0\~ F(t) <F(t) fexp ( -e + ^-^\^A - 1 



+ XN V~ e g| 2 expL-e+ e2(S ^ l \iA ds\ sup F(t + s) 



2 

z(«-«) _ i 

se(0,e) 
We deduce the following contraction estimate, 

li msU p F(t + e) ~ F(t) < -F(t) + (x - 1) (liminf |f*(i + e)| 2 ^) F(t) . 

Hence the one-dimensional Keller-Segel equation (|4.ip is a contraction with rate 1 with respect 
to the Fourier distance d\. D 

15 



Acknowledgments 

Both authors are supported by the bilateral France-Spain project FR2009-0019. JAC is 
partially supported by the projects MTM2008-06349-C03-03 DGI-MCI (Spain) and 2009- 
SGR-345 from AGAUR-Generalitat de Catalunya. Both authors thank CRM (Centre de 
Recerca Matematica) in Barcelona where part of this work was done during the stay of the 
first author as a part of the thematic program in Mathematical Biology in 2009. 

References 

[1] L. A. Ambrosio, N. Gigli, G. Savare, Gradient flows in metric spaces and in the space of 
probability measures, Lectures in Mathematics, Birkhauser, 2005. 

[2] F. Barthe, Inegalites de Brascamp-Lieb et convexite, C. R. Math. Acad. Sci. Paris 324 
(1997), no. 8, 885-888. 

[3] F. Barthe, On a reverse form of the Brascamp-Lieb inequality, Invent. Math. 134 (1998), 
no. 2, 335-361. 

[4] W. Beckner, Sharp Sobolev inequalities on the sphere and the Moser-Trudinger inequality, 
Ann. of Math. 138 (1993), no. 1, 213-242. 

[5] P. Biler, G. Karch, P. Laurengot, T. Nadzieja, The Sir-problem for radially symmetric 
solutions of a chemotaxis model in the plane, Math. Methods Appl. Sci. 29 (2006), no. 
13, 1563-1583. 

[6] P. Biler, T. Nadzieja, Existence and nonexistence of solutions for a model of gravitational 
interaction of particles I, Colloq. Math. 66 (1994), no. 2, 319-334. 

[7] A. Blanchet, V. Calvez, J. A. Carrillo, Convergence of the mass-transport steepest descent 
scheme for the subcritical Patlak-Keller-Segel model, SIAM J. Numer. Anal. 46 (2008), 
no. 2, 691-721. 

[8] A. Blanchet, E. Carlen, J. A. Carrillo, Functional inequalities, thick tails and asymptotics 
for the critical mass Patlak-Keller-Segel model, preprint. 

[9] A. Blanchet, J. A. Carrillo, P. Laurencot, Critical mass for a Patlak-Keller-Segel model 
with degenerate diffusion in higher dimensions, Calc. Var. Partial Differential Equations 
35 (2009), no. 2, 133-168. 

[10] A. Blanchet, J. A. Carrillo, N. Masmoudi, Infinite Time Aggregation for the Critical 
Patlak-Keller-Segel model in R 2 , Comm. Pure Appl. Math. 61 (2008), no. 10, 1449-1481. 

[11] A. Blanchet, J. Dolbeault, M. Escobedo, J. Fernandez, Asymptotic behaviour for small 
mass in the two-dimensional parabolic- elliptic Keller-Segel model, J. Math. Anal. Appl. 
361 (2010), no. 2, 533-542. 



16 



[12] A. Blanchet, J. Dolbeault, B. Perthame, Two-dimensional Keller-Segel model: optimal 
critical mass and qualitative properties of the solutions, Electron. J. Differential Equa- 
tions, (2006), no. 44, 32 pp. (electronic). 

[13] Y. Brenier, Polar factorization and monotone rearrangement of vector-valued functions, 
Coram. Pure Appl. Math. 44 (1991), no. 4, 375-417.. 

[14] V. Calvez, J. A. Carrillo, work in preparation. 

[15] V. Calvez, B. Perthame, M. Sharif! tabar, Modified Keller-Segel system and critical mass 
for the log interaction kernel, Stochastic analysis and partial differential equations, 45-62, 
Contemp. Math., 429, Amer. Math. Soc, Providence, RI, 2007 

[16] E. Carlen, J. A. Carrillo, M. Loss, Hardy-Littlewood-Sobolev Inequalities via fast diffusion 
flows, preprint. 

[17] E. Carlen, M. Loss, Competing symmetries, the logarithmic HLS inequality and Onofri's 
inequality on S n , Geom. Funct. Anal. 2 (1992), no. 1, 90-104. 

[18] J. A. Carrillo, G. Toscani, Contractive probability metrics and asymptotic behavior of 
dissipative kinetic equations, Riv. Mat. Univ. Parma 6 (2007), 75-198. 

[19] D. Cordero-Erausquin, B. Nazaret, C. Villani, A mass-transportation approach to sharp 
Sobolev and Gagliardo-Nirenberg inequalities, Adv. Math. 182 (2004), no. 2, 307-332. 

[20] J. Dolbeault, B. Perthame, Optimal critical mass in the two-dimensional Keller-Segel 
model in R 2 , C. R. Math. Acad. Sci. Paris 339 (2004), no. 9, 611-616. 

[21] B. During, D. Matthes, G. Toscani, A Boltzmann-type approach to the formation of 
wealth distribution curves, Riv. Mat. Univ. Parma 1 (2009), 199-261. 

[22] R. Jordan, D. Kinderlehrer, F. Otto, The variational formulation of the Fokker-Planck 
equation, SIAM J. Math. Anal. 29 (1998), no. 1, 1-17. 

[23] N. I. Kavallaris, P. Souplet, Grow-up rate and refined asymptotics for a two-dimensional 
Patlak-Keller-Segel model in a disk, SIAM J. Math. Anal. 40 (2008/09), no. 5, 1852-1881. 

[24] E. F. Keller, L. A. Segel, Initiation of slime mold aggregation viewed as an instability, J. 
Theor. Biol. 26 (1970), no. 3, 399-415. 

[25] E. H. Lieb, M. Loss, Analysis, Vol. 14, Amer. Math. Soc, Providence, 1997. 

[26] R. J. McCann, Existence and uniqueness of monotone measure-preserving maps, Duke 
Math. J. 80 (1995), no. 2, 309-323. 

[27] R. J. McCann, A convexity principle for interacting gases, Adv. Math. 128 (1997), no. 
1, 153-179. 

[28] T. Nagai, Blow-up of radially symmetric solutions to a chemotaxis system, Adv. Math. 
Sci. Appl. 5 (1995), no. 2, 581-601 

17 



[29] F. Otto, The geometry of dissipative evolution equations: the porous medium equation, 
Comm. Partial Differential Equations 26 (2001), no. 1-2, 101-174. 

[30] C. S. Patlak, Random walk with persistence and external bias, Bull. Math. Biophys. 15 
(1953), 311-338. 

[31] C. Villani, Topics in optimal transportation, Graduate Studies in Mathematics Vol. 58, 
Amer. Math. Soc, Providence, 2003. 



18 



